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ABSTRACT 

Ultracompact X-ray binaries (UCXBs) are thought to have relatively simple binary evolution post- 
contact, leading to clear predictions of their luminosity function. We test these predictions by studying 
the long-term behavior of known UCXBs in our Galaxy, principally using data from the All-Sky 
Survey and the Galactic Bulge Scans with RXTE's PCA instrument. Strong luminosity variations 
are common (and well-documented) among persistent UCXBs, which requires an explanation other 
than the disk instability mechanism. We measure the luminosity function of known UCXBs in the 
Milky Way, which extends to lower luminosities than some proposed theoretical luminosity functions of 
UCXBs. The difference between field and globular cluster X-ray luminosity functions in other galaxies 
cannot be explained by an increased fraction of UCXBs in globular clusters. Instead, our measured 
luminosity function suggests that UCXBs only make up a small fraction of the X-ray binaries above 
a few xlO^^ ergs/s in both old field populations and globular clusters. 
Keywords: binaries:X-ray — globular clusters: general — accretion. 



1. INTRODUCTION 

Ultra-compact X-ray binaries (UCXBs) contain a com- 
pact accretor star (in all known cases, a neutron star 
[NS], though black hole UCXBs are possible) and a 
compact donor star, with an orbital period Porb < 80 
minutes. Such short periods require the donors to be 
hydr ogen-deficient, partially or fully degenerate star s 
(e.g. iRappaport et al.l 119821: iDelove fc BildstenI I2003D . 
UCXB systems are more common in globular clusters 
(GCs), likely due to their formation there by close dy- 
namical interaction s (|Verbuntlll987tlDeutsch et al.ll200Cll : 
llvanova et 10120051 ). 

UCXB systems can be roughly divided into persistent 
(over the decades we have been observing them) and 
transient systems. Transient UCXBs spend the majority 
of the time in a quiescent state with little or no accretion, 
punctuated by occasional outbursts (when they become 
quite luminous, and are more easily detected). The out- 
bursts of known transient UCXBs are rather short, which 
is expected, given the small size of their accretion disks. 

Chandra X-ray observations of elliptical galaxies de- 
tect large number s of X-ray binaries i n GCs (e.g. 
Sarazin e t all 120011: lAng elini et all 120011; 'K undu et atl 
2002; Mi nniti et all [200 4: Jordan et al. 2004). The lu- 
minosity function of GC X-ray binaries differs from the 
non-cluster X-ray binary luminosity function, particu- 
larly below an X-ray luminosity L x = 10'^^ ergs/s where 
GCs have fewer X-ray binaries ( Fabbiano et al.l [20071 : 
iKim et all [200l iVoss et al.l 120091: iZhang et al.l 
UCXBs hav e been suggested to dominate bright GC X- 
ray sources ()Bildsten fc Delo ve 2004), and their high fre- 
quency in GCs has been suggested to explain the dif- 
ferent X -ray lu minosity functions inside vs. outside of 
clu sters (iVoss e t al. 2009). This suggestion depends on 
the iFragos et al.. (|2008i) model of the UCXB luminosity 
function, which cuts off below 5x10^^ ergs s~^ (follow- 

^ Now at International Space University, 1 rue Jean- 
Dominique Cassini, 67400 lUkirch-GrafTenstaden, France 
^ Ingenuity New Faculty; heinke@ualberta.ca 



ing current calculations of the He disk instability line), 
to explain the flat lu minosity function of GCs at low Lx- 
llvanova et al.l (|2008f ) found that UCXBs are unlikely to 
constitute the majority of GC X-ray binaries due to the 
short lifetimes of persistent UCXBs; instead, main se- 
quence donors are preferred for that role. An empirical 
measurement of the UCXB luminosity function (along 
with re-consideration of the appropriate Merit) can un- 
cover the true role of UCXBs among the X-ray binaries 
in elliptical galaxy GCs, and is a target of this work. 

In this paper, we compile histograms of the lumi- 
nosities of galactic UCXBs from the most sensitive sur- 
veys availabl e, and calculate emp irical UCXB luminosity 
functions. In lHeinke et al.l ()2013l ). hereafter Paper II, we 
use the luminosities for individual sources found in this 
paper to discuss their disk stability and evolution. 

2. DATA ANALYSIS 

2.1. Sample selection 

As of January 2013, 13 galactic UCXBs have reli- 
able orbital period measurements (some tentative, see 
Table 1). Four other likely UCXBs show strong indi- 
cations of their ultracom pact nature , prin cipally their 
X-ray/optical flux rat io (iBassa e t aP 12006') . lack of H 
lines in their spectra (iNelemans et al. 2004, 2006), and 
burst characteristics ([Gallowav et al.ii201Qi) . Eight ob- 
jects have been suggested to be UCXBs based principally 
on thei r low (<2 % of E ddington) persistent X-ray lumi- 
nosity ()in't Zand et al.llTOO?. ). To study the UCXB X-ray 
luminosity function, we choose to omit the 8 objects sug- 
gested on this basis alone, as this method of identifica- 
tion would clearly bias the derived luminosity function. 
We also note that a counter-example to the in't Zand 
method-an apparently persistent system below 2 % of its 
Eddington lum inosity, not edge-o n, and with a 2.15 hour 
orbital period (|Engel et al.ll2012D —is now known, H1825- 
331 in the globular cluster NGC 6652. Whether this 
system is a tran sient or not is debatable (jVerbunt et al.l 
ll995tlDeutsch et al.„1998i) . but if it is a transient, it has 
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been active as long as many systems thought to be per- 
sistent, so the discrepancy with the in't Zand method 
remains. 

We cannot be certain that other identification methods 
(which are not uniform) do not introduce luminosity bi- 
ases, but at least they are not obviously biased. We list 
the remaining 17 systems in Table 1, separating them 
into transient and persistent systems. Of these 17, 5 ac- 
creting millisecond pulsars are known transients, while 
the rest are (so far) deemed persistent. 

From the candid ate ultracompacts list e d b y 
lin't Zand et"an (|2007[ ) and iNelemans fc Jonked (|2010f ). 
we exclude 4U 1822-00 ( due to its optical mo dulation 
period of 191 minutes, iShahbaz et all [200l . NGC 
6652 A = H 1825-331 (due to its 2.15 hr optical 
period, lEngel et all 120121 ). NG C 6652 B (due to its 
apparently main-sequence donor Hcinkc et 311120011 and 
recent evidence a gainst the proposed 45-minute period 
lEngel et al.ll2012[ ). and the uj Cen qLM XB (due to its 
strong Ha emission. iHaggard et al.|[200^ . We also omit 
4U 1905-1-000, due to a lack of high-cadence sensitive 
monitoring data (it has been quiescent since at least 
1992; Juett & Chakrabarty 2005). 

2.2. Data selection & analysis 

We took our data from all-sky monitoring surveys 
wherever feasible. The Monitor of All Sky X-ray Image 
(MAX I) detector (Matsuoka ct al. 2009; Sugizaki ct ^ 
120 111) aboard the International Space Station (ISS) pro- 
vides high-sensitivity lightcurves for bright X-ray bi- 
naries anywhere in the sky since August 2009, while 
the Rossi X-ray Timing Explorer's (RXTE) Propor- 
tional Counter Array (PCA, Jahoda ct al. 1996) scans 
of the Galactic Bulge (Swa nk &: Markwardtl 120010 per- 
mit higher-sensitivity lightcurves of X-ray binaries in the 
Galactic Bulge from 1999 to 2012 (and expanded regions 
for shorter times). 

We obtained daily binned MAXI Ga.s Slit Camera data, 
in units of counts s~^ cm~^ in the 2.0-4.0 keV and 4.0- 
10.0 keV energy ranges, from the MAXI websit£0. We 
chose to omit the 10-20 keV energy range, due princi- 
pally to its poorer signal-to-noise ratio, but also because 
we are focused on comparing to Chandra detections of 
extragalactic X-ray binaries below 10 keV. PCA bulge 
scan data were obtained from C. Markwardt's webpagqj 
in the 2-60 keV band. MAXI data collection began on 
MJD 55058 (2009/8/15), with our data spanning approx- 
imately 2 years and 8 months, while there is over 12 years 
worth of PCA bulge scan data for most sources (the bulge 
scans were extended to cover larger regions later). 

We used MAXI data for persistent sources outside the 
bulge scan region, combining the 2-4 and 4-10 keV data 
and errors. We discarded roughly 30 data points per 
source corresponding to dates the shuttle was docked 
at the ISS, and high data points that were clearly at- 
tributable to Sun glints (identifiable as a bright, elon- 
gated source moving over the MAXI map around the 
source over a few days). We also identified data with 
large errors as noisy; the value of the error cutoff was 
determined on an individual basis for each source. This 
amounted to ^0.5-10% of the data from each source. 



Combined MAXI datapoints less than 3 a above zero 
were considered upper limits; we create 3 a upper limit 
points by replacing the datapoint with an upper limit 
three times the error above zero. For PCA data, we used 
a 4 CT detection limit, since the PCA is more sensitive, 
and the crowded bulge leads to increased systematic er- 
rors. For systems known to be transients (e.g. observa- 
tions with sensitive X-ray detectors found them in qui- 
escence) , upper limits were taken to indicate quiescence, 
and were simply replaced with a value of 0. Persistent 
sources were essentially always detected in PCA data. 

Many of these UCXBs exhibit substantial variability. 
The lightcurves of 4U 1850-087, 4U 1728-34, and 2S 0918- 
549 (Fig. 1) demonstrate variation of up to an order of 
magnitude, on timescales from days t o months. Such 
variability was highlighted recently by iMac carone et al.l 
jgoTo) for 4U 1543-40 (the UCXB in NGC 1851). Similar 
variability was noted from, e.g. the UCXBs 4U 1915- 
05by|Simon| (l2005l) and lA 1246-588 bv lin't ZandeFall 
(|2008l) : see below. 

We converted the MAXI data, into 2.0 - 10.0 keV intrin- 
sic luminosities using the Crab nebula as a standard cal- 
ibration source. In the 2-10 keV energy band, an average 
photon flux of 3.12 cts s^^ cm^^ was found for the Crab 
Nebula. Using the Crab's flux in the 2.0-10.0 keV range 
(2.16 X 10-^ ergs cm^^ s^^), we find 6.9 x 10"^ ergs (2-10 
keV)/ct. We assume that our sources have roughly Crab- 
like spectra (power- law of photon index 2.1) within this 
band. Altering the photon index by 1 in either direction, 
or increasing Nh up to ten ti mes larger, alters the unab- 
sorbed fiux by at most 30% (jin't Zand et al.ll2007[ ). For 
the PCA data, the Portable Interactive Multi-mission 
Simulator (PIMMS^ was used to convert photon fluxes 
(counts per s per 5 PCU) into unabsorbed energy fluxes 
in the 2.0 - 10.0 keV range assuming a photon index of 
2 and the Galactic Nh values in Table 1. 

We compute 2-10 keV intrinsic luminosities using the 
best-estimate distances in Table 1. Many of these dis- 
tances come from photospheric radius expansion bursts, 
or from estimates of globular cluster distances, and are 
thus well-constrained. Other distances-for 4U 1543-624, 
4U 1626-67 and four of the five transient systems-are 
relatively poorly constrained (though the location of sev- 
eral of the transients near the Galactic center suggests a 
location in the Galactic bulge). The only distance esti- 
mate for 4U 1543-624 comes from assuming that its M 
is powered by gravitational radiation, and that it follows 
standard white dwarf UCXB evolution. 



3 http://maxi.riken.jp/top/ 

^ http:/ /Iheawww. gsfc.nasa.gov/uscrs/craigm/galscan/main. html ^ http:/ /asc. harvard.edu/toolkit/pimms.jsp 
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Figure 1. MAXI and PCA lightcurves of some persistent UCXB 
systems showing strong flaring; 2S 0918-549 (MAXI), 4U 1728-34 
(PCA), and 4U 1850-087 (PCA). 
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Table 1 
UCXB Information 



Source Location Distance Period Nh MAXI points PCA points 

(kpc) (mins) (10^^ cm"^) Ncood NLimits Ncood ^Limits 



Persistent systems 



4U 1728-34 


Field 


5.2 ±0.8" 


ICSV" 


22.9" 






803 





4U 1820-303 


GC 


7.9 ±0.4' 


11" 


1.6"* 


600 





2013 





4U 0513-40 


GO 


12.1 ±0.6'' 


17= 


0.26'' 


144 


507 






2S 0918-549 


Field 


5.4 ± 0.8^' 


17.4S 


3.0'' 


300 


519 






4U 1543-624 


Field 


7.0' 


18.2' 


3.5'' 


667 


50 






4U 1850-087 


GC 


6.9 ± 0.3' 


20.6^ 


3.9^* 


106 


645 


195 





M15 X-2 


GC 


10.4 ±0.5' 


22.6* 


0.67' 










4U 1626-67 


Field 


g+5 rn 


42m 


1.4" 


739 


108 






4U 1916-053 


Field 


9.3 ±1.4° 


50^ 


3.29 


332 


453 






4U 0614-f091 


Field 


3.2 ±0.5'' 


51?" 


3.0' 


567 


50 






lA 1246-588 


Field 


4.3 ±0.6" 


7 


2.5" 


122 


389 






4U 1812-12 


Field 


4.6 ±0.7" 


7 


15.0"" 






405 





Transient systems 


XTE_J1807-294 


Bulge 


0+4 X 
o_3 3 


40.1^ 


5.6^ 


105 


602 


35 


775 


XTE_J1751-305 


Bulge 


O + 0.5 a 
^-1.3 


42/3 


9.8T 










XTE J0929-314 


Field 




43.6"* 


0.76"= 










Swift J1756. 9-2508 


Bulge 


8±4« 


54.7^ 


54^ 






10 


727 


NGC 6440 X-2 


GC 


8.5±0.4' 


57.3'' 


5.9' 






36 


785 



Note. — Known and suspected UCXBs in our sample, with 
best estimates of their distance, period and A'^^f, plus the number 
of usable data points from MAXI and the PCA bulge scans. Loca- 
tion in the Galactic field, (direction of the) bulge, or in a Globular 
Cluster (GC) is also specified. Errors are ranges fro m indirect es- 
timates; 15% errors on bursts IIKuulkers et al.ll200g) : 5% errors on 
GC distances. Periods supported by only weak evidence have '?'s. 
References: "Galloway ct aL| |[20TO|) : TO arris' ("2010); '^Stel la et al.l 
{1983); ISidb li ct al. (2001); IZurek et al. (2009); -^in't Za nd et al.1 
I l2005al 1: 5|7hong & Wang (201T); ITuett fe'CSiakrabartvl 1 12003) ; 
IWang fc Chakr abarty (2004 ) , distance estimate assumes M driven 
by gravitational radiation; Homer et al.l ( 119961) : * IDieball et al.l 
(|2TO5l); IWhite & Angelinil (|200ll); 'IChaEabart^ ([1998 
lErauss ct al. (2007); ^Yoshida (1993); ? Walter et al. (1982 
• iChurch et al. (1998); 'Brandt ct al. (1992); "Shahbaz ct al] 



( l200l): tP iraino ct al. (1999); "in't Zand et al. (2008 
'I Bassa et a . (2006); "Tarana ct al.. 120061: ^Gallo way (200% 
WarkwardT ct al. ( 20031): IFalanea et al.l (120051') : "TPapitto et S\ 
([2003): /^Markwardt et al.l ( [2003) ; IMiller et al.. (2003) ; 
iGallowav ct al. (2002); ' Jue tt et all (12003) ; HKrimm et"an 
(12007.1 : ' tAltamirano et al.. (.201011 . 
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For four faint transients (XTE J1751-305, NGC 
6440 X-2, XTE J0929-314, and Swift J1756.9-2508) 
and M15 X-2, neither MAXI nor PGA bulge scan 
lightcurves provide a reliable history. For XTE J1751- 
305, we create daily 2-10 keV flux lightcurves of the 
2002, 2005, 2007 and 2009 outbursts from the lit- 
erature jM arkwardt ct al. 2002; Grebc nev et all 120051: 
Swank et al. 2005; Fala nga et al.M2007; ,Markwardt et al] 
20071; [Riggio ct al. 2011), and take other PGA bulge scan 
measurements as evidence of quiescence. For NGG 6440 
X-2, we use the PGA bulge scan obser vatio ns as for other 
sources (as it was det ected in several, IHeinke et aLilSQlOl : 
iPatruno et al.|[2010bD . but exclude the three major out- 
bursts, which were all produced by the oth er known tran- 
sient i n that cluster, SA X J1748.9-2021 (jPatruno et al.l 
I2010at lAltamirano et afll2008l ). As it is a transient, we 
only use detections, and set upper limit observations to 
indicate quiescence. For X TE J0929-314, we use the out- 
burst lightcurve data from lGallowav et all (|2002f ). It was 
clearly observed by the RXTE All-Sky Monitor (ASMfl 
in its 2002 outburst, but was not otherwise detected by 
RXTE ASM, so we consider the remaining ASM his- 
tory since 1996 to be a record of quiescence. For Swift 
J175 6.9-2508, we u s e the 2007 outburst lightcurve data 
from iKrimm et al.l (|2007f) along with the PGA bulge 
scans, converting the Swift Burst Alert Telescope (BAT) 
15-50 keV fluxes to 2-10 keV fluxes using a typical flux 
ratio of 3 (which aligns the BAT and PGA flux estimates 
when they are simultaneous). The 2007 outburst was 
missed by the RXTE PGA bulge scans due to Sun con- 
straints, while other outbursts during the bulge scan era 
have been detected, so we take bulge scan upper limits 
as evidence of quiescence. 

Ml 5 X-2 cannot be separated from AG 211 (X 
2127-1-12) with any X-ray instrument but Ghandra. 
We use one archival and six reported Chandra obser- 
vations of M15 t o inter pret its long-term lightcurves. 
iWhite fc AngelhH (pOOl 's Chandra HETGS spectral fit 
gives Lx{2-10 keV)= 9 x 10'^^ ergs/s. Applying this 
spectral fit to three Chan dra-BRC observa tions in 2001 
( Hannikainen et 311120051 ) and one in 2007 (|Heinke et al.l 
[20091) gives Lx = 1.0 ± 0.1 x 10^^ ergs/s in each case. 
We extract M15 X-2's readout streak spectrum (using the 
standard analysis recipe Q) from archival Chandra AGIS- 
S ObsID 11029 find ing Lx = 7 x 10^^ ergs/s. Finally, 
iSivakoff et all ()20lH) report that M15 X-2 is responsi- 
ble for a major X-ray brightening of M15 in 2011, giving 
Lx = 1.1 X io^"^ ergs/s. We thus suggest that M15 X-2 is 
responsible for similar previous large flares in M15 seen 
by the RXTE ASM. M15 X-2 appears to be usually in 
the range of 7 x 10^^ < < 1 x 10^^ ergs/s, except 
for the bright flares seen by RXTE ASM and MAXL AG 
211 i s typically fainter than M15 X-2, but more vari- 
able (|Hannikainen et al.ll2005l) . Thus we take M15 X-2 
to be between 7 x 10^^ < Lx < 1 x 10^^ ergs/s when 
the total flux from M15 is below 2 x lO'^^ ergs/s. We 
use RXTE ASM data to get better long-term statistics 
of the bright fiares (as only one, very bright fiare was de- 
tected with MAXI). We take a 5-day average of RXTE 
ASM datapoints, and select data with countrate > 1.5 



cts/s, and errors < 0.5 cts/s, as "fiaring", finding 68 such 
datapoints and 915 "normal" datapoints. Although the 
resulting luminosity function is overly simplistic, it seems 
to represent the key elements of M15 X-2's behavior rea- 
sonably well. 

2.3. X-ray luminosity functions 

Galibrated data points were arranged into thirty-one 
luminosity bins between and 10'^^ ergs/s, where all lu- 
minosity values less than 10^^ were placed in the low- 
est bin. The errors for each bin were computed us- 
ing Gehrels' u pper limit approximation, 1 -I- \/N + 0.75 
(|Gehrel£lll98 6V and taken to be symmetric. X-ray lumi- 
nosity functions (XLFs) for each source can be found in 
Fig. 2, with upper limits also plotted for Mj4X/ sources. 
For the five transient sources considered, the XLFs in 
Fig. 2 represent periods of outburst only. 

We create empirical XLFs for the UCXB population 
in two ways: taking one observation for each source 
creates a "snapshot" luminosity function of the UGXB 
population (as from a single observation of our Galaxy) 
and compiling 100 observations of each source creates 
a combined luminosity function over time. Both re- 
quire random sampling from their intrinsic luminosity 
functions, incorporating our detections and upper lim- 
its. For PGA-detected systems, we use the detections 
alone (Table 2). For known transients (Table 3), upper 
limits are assumed to indicate the source is in deep qui- 
escence {Lx << lO'^^ ergs/s), as is typically found by 
deep Chandra or XMM observations. This assumption 
appears valid given both the rapid evolution of transient 
outbursts and recent results monitoring Galactic globu- 
lar clusters for low-luminosity transients (Altamirano et 
al. 2013, in prep). 

For persistent sources with upper limit observations, 
we statistically estimate the true fluxes represented by 
upper limit s usin g the maximum likelihood method of 
lAvni et all ()1980D . This method uses an analytic, re- 
cursive function to estimate the 'true' shape represented 
by input data and upper limits, and was designed to be 
used with relatively sparse and binned astronomical data. 
Avni presents a formula appropriate for data with lower 
limits; since in our case, we deal with upper limits, we 
delineate a brief derivation of our t reatment, modified 
from Avni, in Appendix A (see also iFeigelson &: Nelsonl 
119851 ) . We compute maximum likelihood X-ray luminos- 
ity functions for the persistent sources for which we have 
upper limits (Fig. 3, Table 4). 



® http:/ /heasarc. gsfc.nasa.gov/docs/xte/asm_products. html 
http:/ /cxc. harvard.edu/ciao/threads/streakoxtract/ 
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Figure 2. Observed XLFs for all 17 UCXB sources studied in 
this paper. XLFs of sources with both MAXI and PCA data are 
computed using the higher- quality PCA data. Upper limits are 
plotted (dashed red lines) where relevant (only for MAXI data). 
For transient sources, upper limits are taken to indicate quiescence 
and a luminosity of 0, and are not plotted. 
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Figure 2. cont'd. 
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Figure 3. Inferred X-ray luminosity functions for MAXI sources 
as computed from good detections and upper limits us ing a mod- 
ified v ersion of the maximum likelihood method of lAvni et al.l 
l| 19801 ). The Y-axis is the fraction of all observations in that bin. 
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We randomly selected 100 observations from each 
UCXB's XLF. The 1700 data points were also combined 
to create a combined luminosity function of the popu- 
lation. Since the transients are rarely in outburst, they 
contribute little to the total UCXB luminosity function; 
most values drawn from their distribution are '0 '. Fig.[^ 
displa ys the combined luminosity function, with iGehrelsl 
(fl986l) errors plotted. The y-axis is in units of sources- 
including the quiescent points (not shown), the sum of all 
bins comes to 17. Repeated samples yielded similar his- 
tograms. Samples computed using only good detections 
(rather than the Avni-computed XLFs) for the MAXI 
sources showed only minor differences in bins greater 
than the peak, and the peak location did not change. 




-Sampled Data 
-Errors 



^ 
0.01 




1 10 

Luminosity (10^^ ergs/s) 



Figure 4. Results of extracting 100 data points at random from 
each of 17 sources' XLFs to make a com bined luminosity function 
for UCXBs. Errors from GehreQ 111989 ). Counts in each bin were 
divided by 100 to produce the plot, whose ordinate is in units of 
sources (out of the 17 total used for sampling). 

We fitted the slope of the 100 individual samplings of 
the XLF, and of the combined XLF, with a simple power- 
law in Sherpa, using the C-statistic due to the low (or 
zero) numbers of counts in many bins. We computed 
Gehrels errors on the number of counts in each bin, and 
then put the XLF into a differential, dN/dL format for 
fitting, by dividing the number of counts per bin by the 
bin width. We fit 7 x 10^^< Lx <10^^ ergs/s. Less lu- 
minous bins are significantly incomplete, due to the lim- 
itations of all-sky monitors. (We suffer incompleteness 
in identification of UCXBs in all bins, but that incom- 
pleteness is not obviously luminosity-dependent, and is 
difficult to quantify.) The 100 individual samplings pro- 
duced acceptable power-law fits, with indices of a (for 
dN{L)/dL = fcL"") ranging from 0.92 to 2.57, with a 
mean of 1.66 and error of 0.28 (Fig. 5). The combined 
XLF was poorly fit with a power-law (Fig. 6; index 1.65, 



0.5 0.6 0.7 0.8 0.9 1 

Index of power-law fit 

Figure 5. Histogram of power-law indices a fit to 100 diff'erential 
XLF samplings {dN{L)/dL = kL~"). Typical index errors per fit 
are about 0.3. 

probability 8e-26, reduced statistic 9.6), likely due to the 
limited number of sources producing a "lumpy" XLF. 
However, we did obtain the same power-law index as 
the mean of the individual samples (1.65). The "lump" 
above 10"^^ ergs/s, for instance, is produced by 4U 1820- 
30. Thus a more complex fit would not be meaningful. 

A relevant question is whether the XLF of globular 
cluster UCXBs differs from that of field UCXBs. Un- 
fortunately, our statistics are too small to significantly 
test this question, as we have only five globular cluster 
UCXBs (one of them transient). Although the globu- 
lar cluster UCXBs include the highest luminosity system 
(4U 1820-30), the other globular UCXBs are not unusual 
compared to the field systems, or to other suspected ul- 
tracompact globular cluster systems (e.g. in Terzan 2, 
lin't Zand etalll200l . 

2.4. Comparing to RXTE ASM studies of galactic 
UCXBs 

We note that the recent work of Ivan Haaften et aTl 
(|2012f ) on the luminosity functions for galactic UCXBs, 
using individual RXTE ASM measurements, argues for 
rapid variations of all sources and finds luminosity esti- 
mates that differ from ours. Their analysis is an impor- 
tant contribution, but we have two concerns about the 
analysis procedures. First, systematic errors in the ASM 
data are not considered, which may affect both their in- 
dividual dwell measurements (by increasing the number 
of spurious datapoints) and the countrates averaged over 
all dwells. Two of their sources show negative countrates 
averaged over all dwells, evidence for the existence of sys- 
tematic errors. 

Our other concern is that luminosities approaching 
10'^* ergs/s for 10~* of the observations suggest ther- 
monuclear X-ray bursts, including normal short helium- 
burning bursts or, for lower mass accretion rates, 
intermediate-duration (up to half-hour) X-ray bursts 
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Bin Max 
(ergs/s) 


Fraction of Detections by Source 


4U 1728-34 


4U 1812-12 


4U 1820-303 


4U 1850-087 


M15 X-2 


2.15x10^5 


_ 


_ 


_ 


0.005 


_ 


3.16x10^5 


- 


- 


- 


0.015 


- 


4.64x10^5 


0.001 


- 


- 


0.144 


- 


6.81x10^^ 


0.001 


0.086 


- 


0.395 


- 


1.00x10^'' 


0.001 


0.889 


- 


0.308 


0.931 


1.47x10^'' 


0.001 


0.022 


- 


0.056 


- 


2.15x10^'' 


0.017 


0.003 


- 


0.031 


- 


3.16x10^^ 


0.047 


- 


- 


0.041 


- 


4.64x1036 


0.182 


- 


- 


0.005 


0.025 


6.81x10^'' 


0.347 








0.025 


l.OOxlO^^ 


0.224 




0.001 




0.014 


1.47x10^'' 


0.142 




0.029 




0.004 


2.15x1037 


0.032 




0.034 






3.16x10^7 


0.003 




0.146 






4.64x10^7 






0.367 






6.81x10^7 






0.392 






l.OOxlO^S 






0.033 







Table 2 

X-ray luminosity functions for the four PCA-detected sources and M15-X2, as used in random sampling calculations. All data for these 

sources gave secure detections. 



Bin Max 

(ergs/s) 


Fraction of Detections by Source 


NCG 6440 X-2 


Swift J1756. 9-2508 


XTE J1751-305 


XTE J1807-294 


XTE J0929-314 


1.00x10^5 


0.956 


0.969 


0.996 


0.964 


0.987 


1.47x10^5 










0.000 


2.15x10^5 




0.001 








3.16x10^5 


0.002 




0.000 






4.64x10^5 


0.002 


0.001 




0.001 


0.000 


6.81x1035 


0.026 


0.004 


0.000 


0.011 


0.001 


1.00xl03<5 


0.001 


0.004 


0.000 


0.001 


0.001 


1.47xl03'5 


0.006 


0.001 




0.007 


0.001 


2.15x1036 


0.005 


0.004 


0.001 




0.004 


3.16x1036 


0.001 


0.008 


0.001 


0.001 


0.003 


4.64x1036 




0.007 


0.001 


0.003 


0.003 


6.81xl036 






0.001 


0.003 


0.001 


1.00x103'' 






0.001 







Table 3 

X-ray luminosity functions for the transient sources used in random sampling calculations. All upper limits are taken to indicate 
quiescent observations. Data for sources except for XTE J 0929-314 are from the RXTE PCA instrument; data for XTE J0929-314 are 

from lGallowav et al.l pOOl V 



Bin Max 
(ergs/s) 



lA 1246-588 



1.00x1035 
1.47x1035 
2.15x1035 
3.16x1035 
4.64x1035 
6.81x1035 
1.00x1036 
1.47x1036 
2.15x1036 
3.16x1036 
4.64x1036 
6.81x1036 
1.00x103'' 
1.47x103'' 
2.15x103'' 



Fraction of Detections by Source 



2S 0918-549 4U 0513-40 4U 0614+091 4U 1626-67 4U 1543-624 4U 1916-053 



0.493 
0.283 
0.167 
0.050 
0.008 



0.202 
0.262 
0.226 
0.182 
0.098 
0.031 



0.545 
0.161 
0.175 
0.105 
0.010 
0.003 



0.027 
0.209 
0.476 
0.223 
0.061 
0.005 



0.007 
0.182 
0.525 
0.278 
0.008 



0.021 
0.087 
0.453 
0.420 
0.020 



0.282 
0.377 
0.200 
0.091 
0.050 
0.001 



Table 4 

Inferred X-ray luminosity functions for the seven sources for which MAXI observ ations were used i n our calculations, as computed using 

the modified maximum likelihood method of lAvni et all ITMCl'l . 
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Figure 6. The combined 'population' differential XLF, compiled 
using 100 randomly-drawn data points from each of the 17 sources' 
XLFs. dN/dL was computed for each bin by dividing the effective 
number of sources per bin by the bin width. A power-law of index 
a=1.65 {dN{L)/dL = kL~°') fits the general slope, although the 
XLF has lumps due to individual sources. 

(|Gallowav et al.l[200l liii't Zand et aLl |2005bf). We can 
roughly estimate the effect on 90 second dwells of short 
bursts if we know the recurrence time and a (the ra- 
tio of persistent to burst fluence). 4U 1728-34, for in- 
stance, shows short bursts with decay timescales of 6.3 s, 
recurrence times of 2.5-5 hours, and a averaging 150 
(from 90-300). Thus, each burst will affect one ASM 
datapoint, they will affect 5 x 10~^-10~^ of the dwells, 
and they will increase the flux by typically a factor of 
2. Fainter systems will typically show rarer bursts, with 
larger fractional increases. Intermediate-duration bursts 
and superbursts last longer than an ASM dwell, reach the 
Eddington limit, and are much rarer (e.g. intermediate- 
dura tion bursts appear o nce per 85 days for SLX 1737- 
282, iFalanga et al.l 120081 ). Bursts from UCXBs show a 
range of fluences between these ext remes, which may 
substa ntially affect the observations of lvan Haaften et all 
((Ml) . Clearly more detailed analysis of the unique 
RXTE/ASM dataset is suggested, and may provide deep 
insight into UCXB burst properties. 

3. DISCUSSION 

3.1. Prevalence of strong variability in persistent 
UCXBs 

The origin of substantial (up to a factor of 10) vari- 
ability seen in a num ber of persistent UC XB sources 
is not understood. M accarone et al.l (j2010t ) highlighted 
such variations in 4U 0513-40. Similar nonperiodic 
variations of up to a factor of 10 have been identi- 
fiedin several persi stent UCXBs: 4U 19 15-05 (jSimonl 
[2005I). lA 1246-58 (linT Zand et al.l 1 20081) . 4U 1626-67 
(jCamero-Arranz et al.l |2010() . M15 X-2 (jsivakoff et al.l 



12 011 1), & lin't Zand et"all (|2007[ )'s (strong candidate) 
UCXBs SAX J17 12.6-3739, SLX 1735- 269, & IRXS 
J170854.4-321857 (lin't Zand et al.ll2005b() . We now add 
4U 1850-087, 2S 0918-549, and 4U 1728-34 to the hst. 
Slightly smaller flaring (of up to a factor of 5) is seen 
during times when the per sistent flux is at pa rticu- 
lar values for 4U 061 4-^09 ()Kuulkers et al.l [2010l ) and 
lin't Zand et"aLl (l2007D's strong UCXB candidate 4U 
1722-30. lin't Zand et al.l (|2007 ) discusses an explana- 
tion of the rapid flaring behavior as a tidal instability 
due to a precessing accretion disk (Whitchurst 1988), 
with a possible dependence on mass transfer rate to ex- 
plain its presence/absence in different epochs for, e.g., 
4U 1722-30. Such variations may not exist for high- 
luminosity {Lx > 3 X 10"^^ ergs/s) ultracompacts, with 
their presumably more mas sive donors in shorter or- 
bits ( Bildsten fc Delovell2004 ). though 4U 1820-30's well- 
known cycles illustrate the existence of strong variability 
in some s uch systems. The p revale nce of strong variation 
reinforces iMaccarone et al.l (|2010[ )'s point about the im- 
portance of considering variation among persistent sys- 
tems when identifying extragalactic transients, and thus 
the importance of empirical observations of the luminosi- 
ties of systems of known type for validating luminosity 
function modeling. 

3.2. Luminosity functions of UCXBs 

We have produced an empirical X-ray luminosity func- 
tion from known UCXBs in the Milky Way. The list of 
Galactic UCXBs is certainly incomplete and may suf- 
fer biases (for instance, we are probably missing per- 
sistent and transient sources below 10^^ ergs/s, so we 
advise readers to treat our luminosity function for bins 
below this simply as upper limits). However, it is impor- 
tant to test theoretical calculations against the best data 
available. We have few sources, and typically less than 
one above Lx = 10^'' ergs/s (see Fig. S]), which makes 
comparison of this UCXB luminosity function with those 
of LMXBs in distant galaxies difficult. However, a few 
nearby galaxies have been studied sufficiently to make a 
useful comparison. 

First, we note that the slope of our X-ray luminos- 
ity function (1.66±0.28 for the samplings, 1.65 for the 
cumulative) down to ^ lO'^^ ergs/s is in good agree- 
ment with the predicted slope for ultracompact X-ray 
binaries ( 1.7-1.8, depe nding on the white dwarf entropy) 
found bv lBildsten fc Delovel (|2004l) . This is not a trivial 
result, since the theoretical calculation did not include 
strong variability among persistent systems, transients, 
or a mixture of UCXBs with different formation mech- 
anisms and disk instability lines (needed to explain, for 
instance, the Porb=40-60 minute high-mass-transfer sys- 
tems; see Paper II). 

Our result does not agree with simula tions of the 
UCXB XLF presented in lFragos et"!!! (|2008( ). which have 
similar slopes but cut off at 5x10^^ ergs/s for both per- 
sistent and transient UCXB systems. Our results in- 
dicate that the UCXB XLF continues down to at least 
1 X lO^s ergs/s. The cutoff in lFragos et al.l (|2008l ) is phys- 
ically motivated by the criterion for disk instability of 
irradia ted pure helium disks, calculated in iLasota et al.l 
(|2008D . and assumed to be appropriate for a population 
of UCXBs. Here, we show empirical evidence against 
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this cutoff; in Paper II we identify the physical reasons 
behind the lowered disk instability limit. 

The agreement of the theoretical UCXB slope with the 
XLF slope of several known elliptical galaxies suggested 
that ultracompact binaries dominate the XLF in globu- 
lar clusters (Bildsten &: Deloye 2004 ). However, deeper 
observations have shown that the XLFs of old popula- 
tions break to shallower slopes at luminosities below 10'^'' 
ergs/s, and that globular clusters have a significantly 
shallower slope at low luminosities. These results have 
been shown for the bulge and globular cluster s ystem 
of M 31 (iKong et al.l 120031: iTrudolvubov fc P ricdhorskv 
[200l IVoss fc Gilfanovl 120071) Cen A (IVoss et al.. .20091 



NGC 3379 (|Fabbiano et al.l I2007D. and for comb ined 
studies of several nearby galaxies ([Zhang et al.ll2011[ ) . In 
particular, the similarity of the slope and 10^^ ergs/s 
cutoff of the observed Cen A globular cluster XLF to 
the theoretical UC XB XLF (c ut off at 5 x lO'^^ ergs/s in 
iFragos et aLll2008| ) motivated IVoss et al.l ()2009f ) to sug- 
gest that the Cen A globular clusters were dominated 
by UCXBs. Although this was an excellent idea (given 
the enhanced population of UCXBs in Milky Way glob- 
ular clusters), our UCXB XLF indicates that this is not 
feasible. Other suggested explanatio ns of the difference 
are significantly more complicated (jFragos et al.l 120081 : 
IKim et al.ll2009i : IZhang et al.ll20TTh . 

The steep XLF we find down to Lx = 10'^^ ergs/s, 
combined with the theoretical expectation that this XLF 
will continue to higher Lx and the relatively flat XLFs 
of both bulge and globular cluster XLFs between 10'^^ < 
Lx < 10^'' ergs/s, indicate that UCXBs make up only a 
minor portion of bulge and globular cluster XLFs above 
a few xlO^^ ergs/s. 

4. CONCLUSIONS 

We have constructed luminosity functions for clearly 
identified UCXBs from the best available long-term 
data, using, in order of priority, RXTE PCA bulge 
scans, MAXI lightcurves, RXTE ASM lightcurves sup- 
plemented with Chandra observations, and literature re- 
ports of the luminosities of systems during transient out- 
bursts. Variability by up to a factor of 10 is quite 
common among persistent UCXBs, which is important 
to consi der when identifying tra nsient behavior in other 
galaxies ()Maccarone et al]|2010[ ). We have taken account 
of upper limits in constructin g these luminosity func- 
tions, using the method of Av ni et al.l (jl980) to calcu- 
late a best estimate of the true luminosity function. We 
combined this information to make "snapshot" luminos- 
ity functions (one sample from each source) and a com- 
bined luminosity function using 100 samples from each 
source. Both methods find consistency with a power-law 
of slope a = 1.66 ± 0.28, extending down to 10'^^ ergs/s 
(below which we are incomplete) . 

Our empirical UCXB luminosity function extends 
down in luminosity without a break to ^ 10^^ ergs/s, 
which disagrees with current theoretical luminosity func- 
tions that cut off at 5 X 10"^® ergs/s. The slope, however, 
is in agreement with theoretical predictions from UCXB 
evolution. Comparing our UCXB luminosity function 
to current measurements of luminosity functions in old 
stellar populations and globular clusters in other galaxies 
suggests that UCXBs make up only a small fraction of 
the X-ray binary population in either case. 



This research has made use of MAXI data provided 
by RIKEN, JAXA and the MAXI team, PCA bulge 
scan data provided by C. Markwardt and the RXTE 
team, RXTE ASM data products provided by the RXTE 
team, Chandra archival data, and the ADS. We acknowl- 
edge financial support from NSERC (Discovery Grants to 
COH, NI, and GRS, and an NSERC USRA, Julie Payette 
NSERC Research Scholarship, and Andre Hamer Post- 
graduate Prize supporting MCE), an Alberta Ingenuity 
New Faculty award to COH, a Canada Research Chair 
supporting NI, and the Avadh Bhatia Fellowship sup- 
porting JCG. 

REFERENCES 

Altamirano, D., Casella, P., Patruno, A., Wijnands, R., & van der 

Klis, M. 2008, ApJ, 674, L45 
Altamirano, D., Patruno, A., Markwardt, C. B., et al. 2010, ApJ, 

712, L58 

Angolini, L., Loewenstein, M., & Mushotzky, R. F. 2001, ApJ, 
557, L35 

Avni, Y., Soltan, A., Tananbaum, H., & Zamorani, G. 1980, ApJ, 
238, 800 

Bassa, C. C, Jonker, P. C, in't Zand, J. J. M., & Verbunt, F. 

2006, A&A, 446, L17 
Bildsten, L., & Deloye, C. J. 2004, ApJ, 607, L119 
Brandt, S., Castro-Tirado, A. J., Lund, N., Dremin, V., Lapsliov, 

I., & Syunyaev, R. 1992, A&A, 262, L15-I- 
Camero-Arranz, A., Finger, M. H., Ikhsanov, N. R., 

Wilson-Hodge, C. A., & Beklen, E. 2010, ApJ, 708, 1500 
Chakrabarty, D. 1998, ApJ, 492, 342 
Church, M. J., Parmar, A. N., Balucinska-Church, M., 

Oosterbroek, T., dal Fiume, D., & Orlandini, M. 1998, A&A, 

338, 556 

Deloye, C. J., & Bildsten, L. 2003, ApJ, 598, 1217 
Deutsch, E. W., Margon, B., & Anderson, S. F. 1998, AJ, 116, 
1301 

Deutsch, E. W., Margon, B., & Anderson, S. F. 2000, ApJ, 530, 
L21 

Dieball, A., Knigge, C, Zurek, D. R., Shara, M. M., Long, K. S., 
Charles, P. A., Hannikainen, D. C, & van Zyl, L. 2005, ApJ, 
634, L105 

Engel, M. C, Heinke, C. O., Sivakoff, G. R., Elshamouty, K. C, 

& Edmonds, P. D. 2012, ApJ, 747, 119 
Fabbiano, C, et al. 2007, arXiv:0710.5126 

Falanga, M., Bonnet-Bidaud, J. M., Poutanen, J., & et al. 2005, 
A&A, 436, 647 

Falanga, M., Chenevez, J., Cumming, A., Kuulkcrs, E., Trap, C, 

& Goldwurm, A. 2008, A&A, 484, 43 
Falanga, M., Soldi, S., Shaw, S., et al. 2007, The Astronomer's 

Telegram, 1046, 1 
Feigelson, E. D., & Nelson, P. I. 1985, ApJ, 293, 192 
Fragos, T., Kalogera, V., Belczynski, K., & et al. 2008, ApJ, 683, 

346 

Galloway, D. K. 2006, in American Institute of Physics 
Conference Series, Vol. 840, The Transient Milky Way: A 
Perspective for MIRAX, ed. F. D'Amico, J. Braga, & 
R. E. Rothschild, 50-54 

Galloway, D. K., Chakrabarty, D., Morgan, E. H., & Romillard, 
R. A. 2002, ApJ, 576, L137 

Galloway, D. K., Muno, M. P., Hartman, J. M., Psaltis, D., & 
Chakrabarty, D. 2008, ApJS, 179, 360 

Galloway, D. K., Yao, Y., Marshall, H., Misanovic, Z., & 
Weinberg, N. 2010, ApJ, 724, 417 

Gehrels, N. 1986, ApJ, 303, 336 

Grebenev, S. A., Molkov, S. V., & Sunyaev, R. A. 2005, The 

Astronomer's Telegram, 446, 1 
Haggard, D., Cool, A. M., Anderson, J., Edmonds, P. D., 

Callanan, P. J., Heinke, C. O., Grindlay, J. E., & Bailyn, C. D. 

2004, ApJ, 613, 512 
Hannikainen, D. C, Charles, P. A., van Zyl, L., Kong, A. K. H., 

Homer, L., Hakala, P., Naylor, T., & Davies, M. B. 2005, 

MNRAS, 357, 325 
Harris, W. E. 2010, arXiv:1012.3224 



16 



Cartwright et al. 



Hcinkc, C. O., ct al. 2010, ApJ, 714, 894 

Hcinkc, C. O., Cohn, H. N., & Lugger, P. M. 2009, ApJ, 692, 584 
Heinke, C. O., Edmonds, P. D., & Grindlay, J. E. 2001, ApJ, 562, 
363 

Heinke, C. O., Engel, M. C, Ivanova, N., Pavlovskii, K., Sivakoff, 

G. R., & Gladstone, J. C. 2013, ApJ, submitted (Paper II) 
Homer, L., Charles, P. A., Naylor, T., van Paradijs, J., Auriere, 

M., & Koch-Miramond, L. 1996, MNRAS, 282, L37 
in't Zand, J. J. M., Bassa, C. G., Jonker, P. G., Keek, L., Verbunt, 

F., Mendez, M., & Markwardt, C. B. 2008, A&A, 485, 183 
in't Zand, J. J. M., Cornelisse, R., & Mendez, M. 2005a, A&A, 

440, 287 

in't Zand, J. J. M., Gumming, A., van der Sluys, M. V., Verbunt, 

P., & Pols, O. R. 2005b, A&A, 441, 675 
in't Zand, J. J. M., Jonker, P. G., & Markwardt, G. B. 2007, 

A&A, 465, 953 
Ivanova, N., Hcinkc, C. O., Rasio, F. A., Belczynski, K., & 

Fregeau, J. M. 2008, MNRAS, 386, 553 
Ivanova, N., Rasio, F. A., Lombardi, Jr., J. G., Dooley, K. L., & 

Proulx, Z. F. 2005, ApJ, 621, L109 
Jahoda, K., Swank, J. H., Giles, A. B., Stark, M. J., Strohmayer, 

T., Zhang, W., & Morgan, E. H. 1996, in SPIE Conf. Ser., Vol. 

2808, SPIE Conf. Ser., ed. O. H. Siegmund & M. A. Gummin, 

59-70 

Jordan, A., Cote, P., Fcrrarcsc, L., & ct al. 2004, ApJ, 613, 279 

Juett, A. M., & Chakrabarty, D. 2003, ApJ, 599, 498 

Juett, A. M., Galloway, D. K., & Chakrabarty, D. 2003, ApJ, 

587, 754 

Juett, A. M., & Chakrabarty, D. 2005, ApJ, 627, 926 
Kim, D., Fabbiano, G., Brassington, N. J., & et al. 2009, ApJ, 
703, 829 

Kong, A. K. H., DiStefano, R., Garcia, M. R., & Greiner, J. 2003, 
ApJ, 585, 298 

Krauss, M. I., Schulz, N. S., Chakrabarty, D., Juett, A. M., & 

Cottam, J. 2007, ApJ, 660, 605 
Krimm, H. A., Markwardt, C. B., Deloye, C. J., et al. 2007, ApJ, 

668, L147 

Kundu, A., Maccarone, T. J., & Zepf, S. E. 2002, ApJ, 574, L5 
Kuulkers, E., den Hartog, P. R,., in't Zand, J. J. M., Verbunt, 

F. W. M., Harris, W. E., & Cocchi, M. 2003, A&A, 399, 663 
Kuulkers, E., in't Zand, J. J. M., Atteia, J.-L., et al. 2010, A&A, 

514, A65+ 

Lasota, J.-R, Dubus, G., & Kruk, K. 2008, A&A, 486, 523 
Maccarone, T. J., Long, K. S., Knigge, C., Dieball, A., & Zurek, 

D. R. 2010, MNRAS, 406, 2087 
Markwardt, C. B., Juda, M., & Swank, J. H. 2003, The 

Astronomer's Telegram, 127, 1 
Markwardt, C. B., Pereira, D., & Swank, J. H. 2007, The 

Astronomer's Telegram, 1051, 1 
Markwardt, C. B., Swank, J. H., Strohmayer, T. E., in 't Zand, 

J. J. M., & Marshall, F. E. 2002, ApJ, 575, L21 
Matsuoka, M., Kawasaki, K., Ueno, S., et al. 2009, PASJ, 61, 999 
Miller, J. M., Wijnands, R., Mendez, M., et al. 2003, ApJ, 583, 

L99 

Minniti, D., Rejkuba, M., Funes, J. G., & Akiyama, S. 2004, ApJ, 
600, 716 



Nelcmans, G., & Jonker, P. G. 2010, New Astron. Rev., 54, 87 
Nelemans, G., Jonker, P. G., Marsh, T. R., & van der Klis, M. 

2004, MNRAS, 348, L7 
Nelemans, G., Jonker, P. G., & Steeghs, D. 2006, MNRAS, 370, 

255 

Papitto, A., Menna, M. T., Burderi, L., di Salvo, T., & Riggio, A. 

2008, MNRAS, 383, 411 
Patruno, A., Altamirano, D., Watts, A., & et al. 2010a, The 

Astronomer's Telegram, 2407, 1 
Patruno, A., Yang, Y., Altamirano, D., et al. 2010b, The 

Astronomer's Telegram, 2672, 1 
Piraino, S., Santangelo, A., Ford, E. C, & Kaaret, P. 1999, A&A, 

349, L77 

Rappaport, S., Joss, P. C, & Webbink, R. F. 1982, ApJ, 254, 616 
Riggio, A., Burderi, L., di Salvo, T., Papitto, A., D'Ai, A., laria, 

R., & Menna, M. T. 2011, A&A, 531, A140 
Sarazin, C. L., Irwin, J. A., & Bregman, J. N. 2001, ApJ, 556, 533 
Shahbaz, T., Watson, C. A., & Hernandez-Peralta, H. 2007, 

MNRAS, 376, 1886 
Shahbaz, T., Watson, C. A., Zurita, C., Villaver, E., & 

Hernandez-Peralta, H. 2008, PASP, 120, 848 
Sidoli, L., Parmar, A. N., Oosterbroek, T., Stella, L., Verbunt, F., 

Masetti, N., & Dal Fiume, D. 2001, A&A, 368, 451 
Sivakoff, G. R., Heinke, C. O., Miller-Jones, J. C. A., Altamirano, 

D., Kuulkers, E., & Morii, M. 2011, The Astronomer's 

Telegram, 3393, 1 
Stella, L., Priedhorsky, W., & White, N. E. 1987, ApJ, 312, L17 
Sugizaki, M., et al. 2011, PASJ, 63, S635 

Swank, J., & Markwardt, K. 2001, in ASP Conf. Ser. 251, New 

Century of X-ray Astronomy, cds. H. Inoue & H. Kunieda (San 

Francisco: ASP), 94 
Swank, J. H., Markwardt, C. B., & Smith, E. A. 2005, The 

Astronomer's Telegram, 449, 1 
Tarana, A., Bazzano, A., Ubcrtini, P., Cocchi, M., Gotz, D., 

Capitanio, F., Bird. A. J., & Fiocchi, M. 2006, A&A, 448, 335 
Trudolyubov, S., & Priedhorsky, W. 2004, ApJ, 616, 821 
Simon, V. 2005, A&A, 436, 263 

van Haaften, L. M., Voss, R., & Nelemans, G. 2012, A&A, 543, 
A121 

Verbunt, F., Bunk, W., Hasinger, G., & Johnston, H. M. 1995, 

A&A, 300, 732 
Verbunt, F. 1987, ApJ, 312, L23 
Voss, R., & Gilfanov, M. 2007, MNRAS, 380, 1685 
Voss, R., Gilfanov, M., Sivakoff, G. R., & et al. 2009, ApJ, 701, 

471 

Walter, F. M., Mason, K. O., Clarke, J. T., Halpern, J., Grindlay, 

J. E., Bowyer, S., & Henry, J. P. 1982, ApJ, 253, L67 
Wang, Z., & Chakrabarty, D. 2004, ApJ, 616, L139 
White, N. E., & Angelini, L. 2001, ApJ, 561, LlOl 
Whitehurst, R. 1988, MNRAS, 232, 35 

Yoshida, K. 1993, PhD thesis, Thesis, Tokyo University, (1993) 
Zhang, Z., Gilfanov, M., Voss, R., et al. 2011, A&A, 533, A33+ 
Zhong, J., & Wang, Z. 2011, ApJ, 729, 8 

Zurek, D. R., Knigge, C, Maccarone, T. J., Dieball, A., & Long, 
K. S. 2009, ApJ, 699, 1113 



5. APPENDIX A 

Consider M luminosity bins indexed by n. If our data consisted only of valid detections, the value of the distribution 
function for the n*'' bin, /„, would be straightforwardly given by 

/„ = ^ (1) 

where N{n) is the number of valid detections in the n*'* bin, and J is the total number of detections in all M bins. 

Since wc arc dealing with upper limits as well as valid detections, wc instead need to define for each of the n bins an 
"effective" number of detections, which combines the valid detections and the most likely distribution of the upper 
limits, 

= ^ (2) 

where Neg{n) is an effective number of detections in the n*'* bin. To calculate this effective number, the nature of 
the most likely distribution of upper limits needs to be considered. 
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An upper limit detection with luminosity value L that falls into the fc*'' bin could be random noise or a true detection 
with a luminosity value anywhere from up to the value L. Thus, the true value could fall in any of n bins, with 
n < k. If we consider an individual such bin n, it can contain a portion of the upper limits in all bins with k > n. To 
find the total number of upper limit counts that truly correspond to bin n, T(n), then, we must perform the following 
sum: 



nn) = t^ (3) 

k=n 2^2 = 1 

Where U (k) is the number of upper limit detections which fall into the k*^^ bin. In the maximum likelihood scenario, 
the way that these U (fc) upper limit counts are distributed amongst all lower bins depends on the distribution function 
- that is, the probability bin n would contain a count divided by the total probability for all bins lower than k. We 
must take the outer sum, over all bins with k > n, because upper limit counts from all of these bins contribute to 
"true" counts in bin n. 

Once we know the way in which upper limits are likely to be distributed among the bins, it is a simple matter to 
derive Nas(n): 

NMn) = N{n) + (4) 

k=7l Z^Z = 1 

Where we have simply added the valid detections to the most likely number of upper limits that actually lie in bin 

n. 

Now, returning to equation [5] and substituting equation 2] for Ncg{n), we obtain: 



E^=i^ (5) 



Isolating /„ 



N{n) 



^ -r U(k) 



J _ v^M U(k) 



Finally, we note that we can rewrite the sum in the denominator as follows, using the fact that by definition, all the 
probabilities in the distribution sum to 1: 



k M M M 

z=l z=l z=k+l z=k+l 



z— 

Making the above replacement, we thus obtain 



N{n) 



j_^M U(k) 



This formula is clearly recursive, due to the presence of the distribution function values on the right hand side. 
If wc start with bin Af , the sum in the denominator disappears (which makes sense, as there are no higher bins to 
potentially contribute their upper limits). Once we have Jm, we can use it to calculate Jm-i and so on, proceeding 
from the highest to the lowest bin recursively. 



